On Theory of the Phonon Perturbed Superradiance 
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The paper examines superradiance in impurity crystals in the field of a coherent phonon wave excited by 
two ultrashort laser pulses via Raman scattering processes at the moment of preparation of the initial state of 
an ensemble of emitters. It is shown that by varying the power of the excitation pulses and their mutual di- 
rection of propagation, one can control the superradiance parameters and extract data on the electron-phonon 
coupling constant and its anisotropy. 
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1. INTRODUCTION 



The problem of spontaneous coherent emission (Dicke superradiance in impurity crystals taking into ac- 
count the interaction of emitters, impurity particles, and thermal phonons, has been studied by a number of 
■researchers [2-7]. Weak emitter-phonon coupling generally leads to temperature-dependent decrease in the super- 
radiance intensity [3 El- But with strong adiabatic electron-phonon coupling, for which the energy spectrum of 
the system acquires a set of electron vibrational levels consisting of a zero phonon line and vibrational repetitions, 
.there is a finite probability of emission of superradiance pulses not only in the zcro-phonon transition but also in 
'the vibronic transitions, with the intensity of these pulses, their delay times, and durations strongly depending 
■on the electron-phonon coupling parameter and the temperature of the sample(whose dimensions are smaller 
|than the radiation wavelength) [SJCI. Similar results were obtained in Refs. 3 and 4 using Thomson model 
which describes Brillouin scattering of light in a resonant medium with the participation of a single absorbed (or 
.emitted) phonon. 

Experimentally, superradiance on electron-vibrational transitions of O2 molecules was realized in the polar 
.dielectrics KCL : O2 (see Refs. 9 and 10). Naboikin et al. |Sj studied Dicke superradiance in mixed molecular 
crystals of diphenyl with pyrene on pyrene centers, which typically display weak electron-phonon coupling. 
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2. STATEMENT OF THE PROBLEM 

Obviously, the parameters of superradiance can be effectively controlled by exciting a coherent phonon wave, 
say by a time-dependent Raman process, instead of exciting thermal phonons, which because of their insignificant 
population are ineffective at low temperatures. In this case the experimental setup changes somewhat: to prepare 
the initial superradiance state at time t = there must be, in addition to an ultrashort pulse of frequency wi, a 
second ultrashort pulse of frequency UJ2 propagating at an angle to the first. Then, in addition to optical excitation 
of impurity centers in the sample at frequencies uji and W2, a coherent phonon wave is excited at the difference 
frequency wi — cj2 = and this wave modulates the frequency of the quantum transitions of the emitters (as 
in the case in which only thermal phonons participate) . Changing the direction of propagation of the excitation 
laser pulses and tuning the difference frequency uji — loi to resonance with different phonon modes, one can not 
only control the superradiance dynamics but also use the behavior of the characteristic superradiance parameters 
to determine the features of electron-phonon coupling of localized electron states with the excited phonon modes 
and the anisotropy of this coupling. Note that the role of a coherent phonon wave in the light echo phenomenon 
was studied by Wilson et al. 21] (see also Rcf. 12). 

We study an ensemble of optical emitters - impurity atoms or molecules imbedded in a crystal lattice. A nat- 
ural approach to describing the electron states of impurity particles in resonant phenomena is to use the idealized 
scheme of a two- level quantum system We assume that each such two- level system interacts with radiation 
field and the coherent phonon wave (excited in this case by ultrashort laser pump radiation). The interaction of 
the quantum levels and thermal phonons can be ignored because the phonon modes of the lattice are " frozen" 
if low temperatures are used in the experimental setup. It is legitimate to study a selectively excited coherent 
phonon wave by classical means, i.e., by assuming that the wave has an amplitude q, a frequency fJ, and a phase 
0. 

We write the Hamiltonian of the system of optical emitters interacting with the radiation field and the coherent 
phonon wave as 

H = Hp + Hph + H,r.t + V{t), (1) 

Where 

Hp — loq ^ , 

3 

^ph = y^^fca^(fc)Q(fc), 

k 

H^nt = Mk)a{k){S+ +SJ)+ H.C.], 

v{t) ^xJ2 s-cos{nt + (/)j). 

j 

Here a{k) is the quantized amplitude of the radiation field for the mode (k, s) with the frequency ojk and po- 
larization Gs, 5*^ ' is the pseudospin variable, which obeys the commutation relations for angular momentum and 
describes the jth {j = 1, 2, N) two level atom with energy splitting luq, and gj{k) and A are the electron-photon 
and electron-phonon coupling constants: 



gj{k) — g{k)exp{i'k ■ rj) = —iy^2nc\k\/V{d ■ es)exp{ik ■ Tj), 

With d and the dipole moment and radius vector of the jth particle and V the quantization volume of the 
radiation field, and the phase of the phonon wave at the point occupied by the jth particle. Throughout this 
paper we take h — 1. 
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3. DERIVATION OF THE MASTER EQUATION 



The statistical operator of a system described by the Hamiltonian satisfies the Liouville equation (the 
evolution of the system is studied over time intervals shorter than the time of irreversible dephasing of the 
polarization of the resonant medium) 



which the canonical transformation 



reduces to 



P^p' + U+pU, 

l'^ = [H,nt[t),p'], 



where 



U = exp < i 



J dt'[Hp + Hph + V{t')] 



Hint = ^ gj{k)a{k) <j Sfexp 



i I LJot + —sin(rit + 0j) 



H.c. > exp[—iuj{k)t] + H.c. 



(2) 
(3) 
(4) 

(5) 
(6) 



To identify the contributions of the harmonics whose frequencies are integral multiplies of 17, it is convenient to 
expand the exponential function in ((SJ in a Bessel-function series based on formula 

oo 

^iasinx ^ J„(a)e"^. 

n—~oo 



Then 

Hint ~ 9j{k)Jn j '^(^) ^ [Sj'exp[i{uJo — uj{k) + nVt)t\ x exp{in(j)j) + H.c] + H.c. 

Equation can be transformed in a standard manner into the integro-differential equation 



^ + z[i/™i(i),p'(0)] 



\H,nt{t)AH^nt(f'),(J{t')\\dt' . 



(7) 



(8) 



Finding the trace in the photon- field variables of both sides of Eq. (8), we arrive at equations for the reduced 
density matrix a = Trp/i[p'(i)], which describes the evolution of the quantum states of the impurity subsystem 
only: 

da 



dt 



dTTrph[Hint{t), [H^nt{t -T),p'{t 



(9) 



Here we have used the fact that Trpfi{[Hint{t), p' {0}]) — which follows from the factorization of p' at < = 0: it 
is assumed that initially the system is prepared in a state which density matrix 



p'(0) = a(0)|0)(0|, 



(10) 



where |0)(0| is the density matrix of the radiation field (we may assume that for the frequencies in the optical 
range the radiation field is in the vacuum state |0) with zero temperature). 

Since in the case of superradiance the field has the shortest correlation time (proportional to the reciprocal of 
optical-frequency bandwidth), in a certain sense the field can be interpreted as a reservoir for the atomic system. 
The interaction of this system with such a wide band reservoir rapidly wipes out any memory of the system about 
its past, with the result that the temporal behavior of the optical emitters is Markovian. This justifies replacing 
p'(t — t) by p'{t) in ^ and extending the upper limit of the r-interval to oo. Moreover, we assume that p'{t) 
can be written as 

p'(t)-a(t)|0)(0| + Ap', 
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where Ap' is at least of order Hint- Then to second order in Hint we obtain a closed equation for the reduced 
matrix density: 

% = ~J^ dTTrph[H,nt{t),[H^nt{t~T),a{t)\0){0\]]. (11) 

After substituting the expression Q for Hint into the righthand side of Eq.(ll), expanding commutators, taking 
the trace, and integrating with respect r we arrive at 



— — ^ '^"^ [fl I "^"^ I \ ^^P^"^^"^ ~ m)Qt]exp{in4)j)exp{-im(j)i) x [g*{k)gi{k)[S^ , a] 



da 

j.l,{k},n,m 

iP 



luq — Lu{k) + mCl 

iP 



+ ttSIluq — uj{k) + mfl] 



luq + uj{k) + mCl 



+ TrS[uJo + U!{k) + mfl] 



+H.C. = -A{t)a. (12) 

In the discussion that follows we ignore the principal part (denoted by P), which yield only a small imagi- 
nary term, i.e., causes a small shift in the resonant frequency, and the nonrcsonant terms containing the factor 
S[uJo + uj{k) + mfi]. 

The characteristic time it takes the optical emitters to go into a coherent state is proportional to ||A||^^ (here 
stands for the operator's value in frequency units) and is much longer than the period T = 27r/il of the rapid 
harmonic "jitter" superimposed on the less rapid evolution process of spontaneous coherent emission. Hence it 
is natural to decompose the real motion, described by the density matrix a, into slow motion averaged over the 
period T and rapid "jitter" . It is convenient to realize the decomposition procedure by employing the method of 
"temporal projection" operators 

a{t)^a{t) = ^ dta{t) = P'<j{t), (13) 

where the projection operator P* performs the averaging of rapidly varying quantities over the period T. Next 
we introduce the operator Q* = 1 — P*. Then the decomposition of the density matrix into the slowly and rapidly 
varying parts can be written as 

CT = P*cr + Q*cr. (14) 

Successively applying the operators P*and Q* to Ec^. (12), we arrive at two differential equations instead of 
one [T^: 

^ = -P*A(i)PV - P*A(i)gV (15) 

9ll ^ -Q*A(t)PV - Q*A(i)QV (16) 
at 

On the basis of Eq. (16) we can formally represent the solution for the rapidly varying (oscillating) part of the 
density matrix in the form 

= - / dt'Q'' K{t')d{t') - f dt'Q'' K{t')a{t') (17) 



Using (17) in Eq. (15) for the slowly varying part of the density matrix, we see that the order of smallness in A 
of the P*AQ*cr is no less than two (or no less than four in Hint), so that the term can be ignored. As a result, 
averaging Eq. (12) over the "jitter" reduces, according to (15), to averaging the operator A(t) over the period 

'^^^Y.Jni^ exp[in{<P, - 00] X E l5(^)l'exp[zk • (r, - r,)] x 5{u:^ - uj{k) + nVL)[S-d, S+] + H.c. (18) 
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This expression constitutes the most general form of the master equation for superradiance in the presence 
of a coherent phonon wave. However, to use this equation that can be compared with experimental data, we 
introduce a "coarsening" procedure. The point here is that real samples containing emitters are much larger that 
the radiation wavelength: koR ^ 1 (here fcg = lo/c and R = max\rj — r;|). Hence exp[ik ■ {rj — r/)] must be 
averaged over the ensemble of particles. We apply a similar procedure to the factor exp[in((j)j — 0/)], assuming 
that the phase difference S(j)ji — (f)j ~ (pi is uniformly distributed over the interval to 27r with density l/27r. The 
result of ensemble averaging, 

ea;p[zk.(r, -rz)] =r(|k|). 

Depends on the shape of the sample, while averaging over the phase spread 5(j)ji results in the exponential function 
exp[inS(j)ji] being replaced by the Kronecker symbol: 

1 f^^ 

— / 6(j)jiexp{in5(j)ji) = Sn.o- (19) 

Allowing for (19) and passing in (18) from the sum over k to an integral with respect to the wave vector k and a 
sum over the polarization s of the photon field, 



we obtain 




where 7 = 2\d\^ uj^ /3c^ . 

Thus the dependence of n-phonon processes (n=l,2,...) on the randomly distributed phase shifts S(pji (with 
constant density ^ must lead to averaging their contribution to zero when superradiance is initiated in large 
samples. The entire modulation effect of a coherent phonon wave on optical centers in this case resides in the 
factor jQ{X/il), which renormalizes the electron photon coupling constant. This factor is similar to the Debye- 
Waller factor, which characterizes the reduction in the unshifted component (the zero-phonon line) of an optical 
transition when elastic scattering of thermal phonons is taken into account. However, in contrast to the case of 
thermal phonons, a coherent phonon wave does not lead to homogeneous dephasing of the polarization of the 
optical centers. 

4. RESULTS AND DISCUSSION 

The standard form 14 of the operator part of Eq. (20) makes it possible to immediately write an (approxi- 
mate) expression for the intensity / of superradiance under conditions in which all emitters are initially completely 
inverted and the transverse macroscopic polarization is zero, i.e., Trp[Sj(j{0)] = ^ and Trp[S^ Sj^ (t{0)] = 0: 

where the renormalized delay time and duration of the superradiance pulse are defined as 

to = TRlnN = t'-^^ J^^ 

(here t'^^ = [7r(A;o)A^]^^) and t^"* are the values of the parameters in the absence of a coherent phonon wave). 

The renormalization factor J§{X/Q), which enters into the expression for to, td, and / because of the modu- 
lation of the quantum states of the emitters by a coherent phonon wave, leads (since Jq (A/fi) < 1 for X/il ^ Q) 
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to an increase in the time of formation of a superradiance pulse, the spreading of the pulse, and a decrease in the 
peak intensity. The advantage of using a coherent phonon wave is that not only does it insure effective control 
over the superradiance process but it also allows extracting information about the strength of the coupling of 
the electron states of the emitters and the various phonon modes thanks to the possibility of selectively exciting 
these modes (and simultaneous "freezing" of thermal phonons caused by low temperatures). These facts are 
illustrated by Fig. 1, which depicts the time dependence of the intensity / for the different values of the total 
energy W of the excitation pulses. Here we have allowed for the fact that the dimensionless parameter X/ft can be 




The intensity I of Dicke superradiance (normalized to Nujo/'It'^) as a function of r = t/r^ for three values 
of the energy W of the excitation laser pulses: 1), W = 0; 2), W = 1.25 x 10^^ J; and 3), W = 2.25 x 10"^ J 
(iV= 10'**). 

written as X/fl = qa^j2Sm^l/h^ where go — AW is the phonon amplitude, S* is a constant (for the oig-mode 
naphthalene, = 1385 cm^^, S — 0.01, A — 0.3 x lO^^^dyn^^, and m ~ lO'^mp, with rup the proton mass). 
Figure 2 illustrates the dependence of a superradiance pulse on the electron-phonon coupling constant A/J7 for 
the fixed value of W). Obviously, in the latter case it is possible to study the anisotropy of the electron-phonon 
coupling by varying the direction of propagation of the phonon wave (i.e, by changing the direction in which the 
excitation laser pulses are applied). 
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